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ABSTRACT 

Some binary evolution scenarios to Type la supernovae include long-period binaries that evolve to 
symbiotic supersoft X-ray sources in their late stage of evolution. However, symbiotic stars with steady 
hydrogen burning on the white dwarf's (WD) surface are very rare, and the X-ray characteristics are 
not well known. SMC3 is one such rare example and a key object for understanding the evolution of 
symbiotic stars to Type la supernovae. SMC3 is an eclipsing symbiotic binary, consisting of a massive 
WD and red giant (RG), with an orbital period of 4.5 years in the Small Magellanic Cloud. The 
long-term V light curve variations are reproduced as orbital variations in the irradiated RG, whose 
atmosphere fills its Roche lobe, thus supporting the idea that the RG supplies matter to the WD at 
rates high enough to maintain steady hydrogen burning on the WD. We also present an eclipse model 
in which an X-ray emitting region around the WD is almost totally occulted by the RG swelling over 
the Roche lobe on the trailing side, although it is always partly obscured by a long spiral tail of neutral 
hydrogen surrounding the binary in the orbital plane. 

Subject headings: binaries: eclipsing - binaries: symbiotic - circumstellar matter - stars: individual 
(SMC3) - stars: winds, outflows - white dwarfs 
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1. INTRODUCTION 

SMC 3 is a symbiotic star consisting of a massive 
white dwarf (WD) and an M giant. Its orbital period 
is about 4.5 years. It is a supersoft X-ray source (SSS) 
(RX J0048. 4 -7332) that has been under observation 
for 20 years (IJordan et al. 1 1199ft IMiirset et al. 1 119971; 
Kahabka 1 120041: iKahabka fc Haberl I l2006t lOrio et al. I 
2007t iSturm et al. 1120111 and references therein). The 
WD appears to be very massive because SMC3 is 
the strongest X-ray source among symbiotic binaries 
(jMiirset et al. I 1 19971) . In addition, a massive WD 
(Mwd > 1.18M©) is suggested from the res ults of X-ray 
spect rum fittings with model atmospheres ()Orio et al. I 
|2007|) . The 20 years observation period is quite long com- 
pared with the durations of the supersoft X-ray phases 
of classical novae harboring massive WDs, 0.7 years for 
1.0 M Q , 2.7 ye ars for 0.8 M^, and 8.9 years for 0.6 M 
(see Table 2 in lKato fc Hachisu I [l99l . Th erefore, it is 
very u nlikely that SMC3 is a nova remnant. lOrio et al. I 
( 2007) concluded that SMC3 is a symbiotic binary with a 
steady hydrogen-shell-burning WD because no nova out- 
bursts have been recorded in the last 50 years despite 
continuous optical monitoring of the SMC. 

Steady hydrogen-shell-burning WDs have attracted at- 
tention in relation to Type la supernova (SN la) pro- 
genitors because these WDs accumulate accreted mat- 
ter, causing their masses to increase. In the single de- 
generate (SD) scenario to SNc la in which an accret- 
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ing WD grows in mass to reach the Chandrasekhar 
mass and explodes as an SN la , there are two well- 
known paths of binary evol ution (|Li fc van den Heuvel I 
[1^9^lHachlsu et al. Ifl999allbl ): one is the main-sequence 
(MS) channel, for those with an MS companion, and 
the other is the symbiotic channel, for those with an 
red giant (RG) companion. In these evolutionary chan- 
nels, a typical binary evolves in its final stage from the 
wind phase (steady hydrogen-shell-burning with opti- 
cally thick winds) through the supersoft X-ray phase 
(steady hydrogen-shell-burning with no optically thick 
winds) to the recurren t nova phase (interm ittent hydro- 
gen shell burning, e.g.. lHachisu et aLll2010f ) . In the sym- 
biotic channel, progenitor binaries are expected to evolve 
into symbiotic stars wi th a massive WD w ith steady 
hydrogen-shell-burning (Hachisu et al.l [20121 ) . They are 
potentially strong SSSs. If they are detected, we use 
them as testbeds for the SN la progenitor scenarios. Un- 
fortunately, symbiotic X-ray binaries with a massive WD 
are very rare. SMC3 is one such rare example. 

It has been suggested that the supersoft X-ray fluxes of 
galaxi es are much weaker than expected fo r the SD sce- 
nario jOrio et al. H2010UDi Stefimoll2l)l7)i3) . IDi Stefano I 
(2010b]) emphasized that the double degenerate (DD) 
scenario for SNe la, in which two WDs merge and then 
explode as an SN la, predicts a large number of symbi- 
otic stars before the second common envelope evolution 
that produces a DD system. Therefore, the lack of su- 
persoft X-ray flux is a problem common to both SD and 
DD scenarios. Gilf anov fc Bogdan I ((2010) claimed that 
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weak supersoft X-ray emission in early-typ e galaxies is 
negati ve evidence against the SD scenario. lOrio et al~l 
(2010) pointed out that most SSSs in M31 are tran- 
sient, and steady-burning sources are rare, possibly be- 
cause of low duty cyc les of the SSS phas e s in these 
sources. The wor k bv iGilfanov fc Boedanl (12010ft was 
also criticized bv [Hachisu et all (120101). Katol (1 2 Oil), 
and iWang fc Hani (|2012t ). because Gilfanov & Bogdan 
adopted not only overestimated numbers of expected 
SSSs but also unrealistically large X-ray fluxes for sym- 
biotic SSSs. If we use a reliabl e number of symbiotic 
SSSs based on the SD scenario (jHachisu et al. I Il999a|) 
and realistic X-ray fluxes, we can statistically explain 
the weak supersoft X-ray emission in early-type galax- 
ies, consequently, this strongly su pports the SD mode l 
as a promising scena rio for SNe la (jHachisu et al.ll2010h . 
iDildav et al. I (j2012f ) reported a detailed spectral analy- 
sis of the SN la PTF llkx and concluded that it had a 
symbiotic nova progenitor, mainly because o f a com plex 
circumstellar environment. Chiotelli s et al. I (|2012t ) also 
suggested that Kepler's SN (SN 1604) is Type la, and 
its progenitor was a symbiotic star consisting of a WD 
and an asymptotic giant branch (AGB) star because of 
strong interaction with circumstellar matter. 

However, the nature of symbiotic steady-burning X- 
ray sources are yet to be fully understood; for exam- 
ple, their WD masses, RG companion masses, spectral 
types, orbital periods, wind mass-loss rates from both 
stars, and nebular emissions, as well as their origins and 
evolutionary paths to the present binary state. In sym- 
biotic stars, a hot WD may be embe dded in a RG's cool 
wind that obscures supersoft X-ra ys (jHachisu et alj20lol : 
iKato I |20H iNielsen et aH 120121) . For example, V407 
Cyg, a symbiotic classical nova, showed a very weak su- 
persoft X-ray flux in the later phase of its 2010 outburst 
because the WD is deeply embedded i n a thick enve- 
lope originating from the cool giant wind (jSchwarz et al. I 
l201~lTlNelson et al. Il20ll IHachisu fc Kato II2012D . Inter- 
estingly, strong symbiotic SSSs have been discovered so 
far only in the galactic halo (AG Dr a) or in the Sma ll 
Magellanic Cloud (SMC3 and Lin358: IOrio et al. Il2007l ). 
This suggests that a low-metallicity environment is an 
important clue to the detection of supersoft X-rays from 
symbiotic systems. To examine such a possibility, we 
need to study individual symbiotic SSSs in detail. 

In symbiotic binaries, a cool wind blows from the RG, 
and a hot wind blows from the WD. Radiation at var- 
ious wavelengths may originate from different areas of 
the binary. Researchers have become interested mainly 
in SMC3's X-ray nature, but its binary nature has not 
yet been fully understood. In relation to binary evolution 
scenarios to SNe la, it is very important to know why the 
supersoft X-ray fluxes are so weak, although their WDs 
are believed to be very massive. In this work, we pro- 
pose a binary model that explains both the long-term 
V magnitude variations and observed X-ray eclipses. In 
Section [2l we summarize the observational properties of 
SMC3 and propose our basic model of it. In Section[3l we 
present our model light curves for long-term variations 
in the V magnitude. Section [4] presents our X-ray eclipse 
light curve model. Discussion and conclusions follow in 
Sections [5] and O respectively. 

2. THE MODEL OF SMC3 - A SYMBIOTIC X-RAY SOURCE 




Fig. 1. — Schematic view of the binary in the V band. The RG 
fills its Roche lobe, and its hemisphere is irradiated by the hot WD 
(right: black dot). Distribution of the RG's surface temperature 
is indicated by red color gradient on a linear scale (white-red side 
is the highest temperature: 5400 - 5500 K). The temperature of 
the RG on the non-irradiated side (dark red) is the parameter Trq 
(3100 - 3500 K). We assume a steady hydrogen-shell- burning WD. 
The WD and its surrounding hot gas are invisible in the V band 
because of their high temperature. The RG tail of neutral hydrogen 
is unseen because it is transparent in the V band. 

Morgan (1992) reported low-resolution optical spectra 
of SMC3. The RG is classified as spectral type M0. The 
spectrum is typical for a high-excitation s-type symbi- 
otic star and exhibits Balmer, He I, and He II emission 
lines as well as Raman-scattered O VI 6825 and 7082 
features. The spectra also show an [Fe X] emission line, 
an unusually high ionization line for a s ymbiotic star. 
Optic al spectra were also published by iMurset et al. I 
(1996) for 1994 Dec 20 - 23 [orbital phase of 4> = 0.26, 
see Eq uation (l2l) in Sectio n [3] below for the ephemeris], 
and bv IQrio et aL I (pOOl for 1994 Oct 17 (</> = 0.22). 
IQrio et al. I (|2007D remarked that coronal lines [Fe IX] 
and [Fe X] in the optical spectra indicate a photoioniz- 
ing source at a very high temperature of T ~ 5 x 10 5 
K. One IUE UV spectrum (SW P 47572L) was taken on 
1993 April 30 (JD 2449107.5) (jVogel fc Morgan I fl99l 
(4> = 0.89, i.e., X-ray ingress). This spectrum ap- 
ears typical for a sy mbiotic star of moderate excitation 
Jordan et al. 1119961 ). 

Near-infrared (IR) / magnitude light curves (shown 
later in Figure [3]) taken by the OGLE II and III projects 
show an irregular variation with a p eriod about 110 days 
(jKahabka 112001 iSturm et al. 11201 ll) which is attributed 
to pulsation of the RG bv lKahabka I (|2004D . A blue filter 
light curve was secured by the MACHO project for as 
long as 6 years (also shown later in Figure [3]); it shows 
a long-term variation with a period of ~ 1600 days simi- 
lar to the X-ray variation (jKahabka Il2004t ISturm et al. I 
2011), superposed by a short-term variation similar to 
that of the I magnitude light curve. 

Therefore, as a basic model of SMC3, we assume a 
binary consisting of a massive WD and an irradiated RG 
as plotted in Figure Q] (which is explained in detail in 
Section 12. 1|) . Table Q] gives the assumed masses of the 
WD and RG, the separation a, and the effective radii of 
the Roche lobes for the RG and WD. Here the binary 
separation a is calculated from Kepler's third law, a 3 = 
G(P rb/27r) 2 (M\yD + -Mrg), assuming a circular orbit 
and P rb = 1634 days (as determined in Section [3]). 

2.1. The RG Component 

Near-IR photometry of SMC3 is reported as (J) = 
11.935, (H) = 11.035, and (K 8 ) = 10.8 at JD 2451034.7 
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TABLE 1 

Binary Parameters 81 



WD mass 


RG mass 


a 


R RL (RG) b 


Krl(WD) c 


K lvl @) 


V vl <3) 








1.35 


3.0 


. 953 


428 


298 


1.35 


2.0 


. 873 


361 


302 


1.35 


1.5 


. 828 


321 


306 


1.35 


1.0 


. 776 


274 


314 


1.35 


0.8 


. 753 


252 


320 


1.35 


0.6 


. 729 


227 


329 


1.35 


0.4 


. 703 


197 


342 


1.20 


0.8 


. 735 


253 


305 


1.10 


0.8 


. 723 


254 


294 


1.00 


0.8 


. 710 


255 


283 



Wc assume P m 



1634 days and a circular orbit. 



D Effective Roche lobe radius of the RG (Egslcton 1983) 
c Effective Roche lobe radius of the WD ( EgKlcton 1983) 



(|Phillips 1120071) (hereafter K s = K). This date is close to 
mid-eclipse, = 0.07. These JHK magnitudes represent 
the RG magnitudes well. Then, we obtain reddening- 
corrected magnitudes of Jo = 11-848, Hq = 10.977, and 
K Q = 10.765 with Aj = 0.874 x E{B - V) = 0.087, 
A H = 0.589 x E(B - V) = 0.058, and A K = 0.353 x 
E(B -V) = 0.353 x 0.099 = 0.035. Here we use th e ex- 
tinction coefficients given by iCardelli et al. I (I1989D and 
the extinction E(B-V) = 0.099 given by[Hilditch et al. I 
<f2005h . In our light curve fittings, on the other hand, the 
RG magnitudes are obtained from the light curve min- 
imum; i.e., when the RG is in front of the WD, we see 
the non-irradiated side of the RG. From Figure [3] we get 
/ ~ 13.34 and V ~ 15.35. Thus, the reddening-corrected 
magnitude becomes Iq = 13.193 and Vq — 15.043, with 
A r = 1.48x0.099 = .147 and A v = 3.1 x 0.099 = 0.307 
(jGardelli et al. lfl989h . 

The RG luminosity is estimated using the relationship 
betwe en BC(K) and ( J — K )g obtained for SMC vari- 
ables (|Bessel fc Wood Ill98l . 



BC{K ) : 



:0.60 + 2.65(J- 
for 0.6 < (J 



K) -0.67(J 
-K) < 1.5. 



K)l 



(1) 



We obtain BC(K) = 2.68 for (J - K) = 1.083. The 
bolometric magnitude becomes -M DO i = 10.765 + 2.68 — 
18.91 = -5.46, i.e., L bo i = 12000 L Q . Here we use 
A'^boi(Sun) = 4.7 5 for the absolute bolo metric luminos- 
ity of the Sun. iBessel fc Wood I (|1984[ ) also presented 
another relationship, BC(K) vs. (V — K)q, from which 
we get BC(K) ~ 2.7; this is quite consistent with the 
estimate obtained from Equation (JXJ) . 

The effective temperature of the RG is estimated using 
the color-temperature relatio n of a model atmosphere for 
low-metallicity RGs given by lKucinskas et al. I (|2006l ) to- 
gether with a metal abundance of [M/H]= — 1. We get 
T off ~ 3800 K for (J - K) = 1.083, T cff - 3600 K for 
(V-K)q =4.278, andT cff - 3700 K for (V-I) = 1.85. 
These color temperatures are consistent with each other 
within the accuracy of the model temperature (±100 K). 
They are a lso consis t ent w ith the spectral type M0 es- 
timated by iMorganl (|1992f ) because for Galactic bright 
gia nts, the spectral type MO co rresponds to 3700 - 3800 
K (|Straizvs fc Kuriliene I [l98l . The radius of the RG 
becomes 250 R Q for T eB = 3800 K, 270 R Q for 3700 
K, and 280 R e for 3600 K, with Stefan-Boltzmann's 
law L = 47ri?p G crT c 4 ff . Comparing these values with the 




Fig. 2. — Schematic polar view of SMC3. Blue circular region 
represents a hot ionized plasma centered on the WD that scat- 
ters supersoft X-ray photons from the WD. The RG almost fills 
its Roche lobe and has a large trailing tail of neutral hydrogen 
that absorbs supersoft X-rays. Arrows denote the orbital motion of 
each component. Ionized and neutral regions are plotted by dotted 
and solid lines, respectively. This schematic fi gure was made us- 
ing information from numerical ca lculations (Sawada et al~l ll986l : 
Mohamcd & Podsiadlowski 2012). Contours denote the density 
taken from lSawadlTeTldl""PIli)86l ). Size and shape of each part are 
not important in our X-ray light curve model except for the height 
(thickness) from the orbital plain. 



TABLE 2 

Light Curve Model of X-Ray Eclipse 



Subject 



Model 1 Model 2 



Duration of total X-ray eclipse Arj> 
Duration of total X-ray eclipse (day) 
Delay of X-ray mid-eclipse Atf> 
Delay of X-ray mid-eclipse (day) 
RG radius 

Minimum height of the tail 
Radius of the hot nebula -Rime 



0.22 

363 day 
0.075 
123 day 
245 Rq 
48 R @ 
100 R Q 



0.19 
304 day 
0.075 
123 day 
199 Rq 
31 Rq 
100 Rq 



Roche lobe radii in Table [T] (i.e., 252 R e for a 0.8 M Q 
RG and 274 R Q for a 1.0 M RG), we see that the es- 
timated RG radii are quite consistent with the effective 
Roche lobe radii of the RG. Thus, we may safely assume 
a Roche-lobe-filling RG (Figured]). 

The RG is a semiregular variable with a pulsation pe- 
riod of 110 days. SMC3 is located in the upper part of the 
period - K magn itude diagram for SMC OGLE variables 
(jlta et al. 1120041 ). In this diagram, SMC3 is located at 
the border between sequences A + and B + , both of which 
are classified as "less regular pulsation AGB variables." 
The pulsation mode is probably not fundamental, but 
is likely the first/second or third overtone. This is con- 
sistent with the irregular / band variation with a small 
amplitude of Al ~ 0.05 mag. 

2.2. The Hot Component: Supersoft X-ray Source 
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As briefly introduced in Section [TJ SMC3 (=RX 
J0048.4— 7332) is an exceptional ly bright SSS among 
symbiotic stars (see Table 1 in iMiirset et al. I H997T) . 
and has been observed with ROSAT, Chandr a , and 
XMM-Newton (IJordan et al. I Il996t IKahabkal 120041; 
IQrio et al. I 120(11" iSturm et al. I 120111 . and references 
therein). IJordan et al. "(j l996fl analyzed ROSAT spectra 
combined with a UV spectrum from IUE. They found 
that a very hot Wolf-Rayet-type model atmosphere fit 
the spectra reasonably well, but blackbody and hydro- 
static plane-parallel models did not. In their analysis, 
the wind mass-loss rate is several times 10 _6 M Q yr _1 
and the central stellar surface is hidden by stellar winds. 
However, they combined the ROSAT (<fi = 0.65: out 
of eclipse) and the IUE spectrum (<f> — 0.89: ingress), 
which are not simultaneous, so they may have overesti- 
mated the wind mass-loss rate in order to reproduce the 
smaller UV flux (due to the half-eclipsed state) with only 
the wind absorption. 

The supersoft X-ray count rates for SMC3 decreased 
by a facto r of 50 over mo re than 300 days starting in 
May 1993 (jKahabka II2004D . Kahabka presented a model 
of this X-ray decrease in which the WD orbits within a 
massive wind from the RG. The wind is highly ionized 
but still contains neutral hydrogen. The column den- 
sity of neutral hydrogen (iVn) in the line of sight changes 
periodically with the orbital motion of the WD. As a 
result, the X-ray count rate decreas es when the WD is 
on far side as seen from the Earth. IQrio et al. I ()2007f ) 
criticized this scenario because this model predicts neb- 
ular emission lines at the high mass-loss rate of the RG 
wind they assumed, but no such lines were observed in 
the optical spectra. Orio et al. argued that the values 
of TVh appeared to be unchanged or only slightly higher 
in the 2003 eclipse, which seemed to be consistent with 
a real eclipse caused by the RG itself, rather than with a 
significant contribution by additional obscuration caused 
by an RG wind. 

Recently, ISturm et al. I (|201lD analyze spectra newly 
obtained with XMM-Newton EPIC-pn, using a black- 
body model, in both high and low X-ray flux stages. 
Their spectral fits do not yield good results for a model 
of changing neutral-hydrogen column density because the 
spectra in low X-ray flux stages cannot be reproduced. 
Therefore, Kahabka's model was not supported. Instead, 
they suggested Compton scattering by fully ionized mat- 
ter that is distributed azimuthally unevenly in the orbital 
plane. This requires the RG cool wind to completely ion- 
ize up to the RG surface because the neutral-hydrogen 
column-density does not increase. However, such a situ- 
ation is v ery unlikely bec a use Raman-scattered light was 
observed (jMorgan Ill992t ). ISturm et al. I (|2011[ ) also sug- 
gested a temperature variation model for the eclipse, but 
this seems to be inconsistent with the light curves. If the 
temperature varies while the total luminosity remains 
constant, the V magnitude light curve should be anti- 
correlated with the X-ray variation, but there is no indi- 
cation of such behavior. Although these authors misun- 
derstood and rejected, the rest of their models, i.e., the 
varying emitting area model with constant temperature 
an d column d e nsity, corresponds to the "real eclipse," 
as IQrio et al. I (|2007l ) already suggested, and reproduces 
very well the spectra at both high and low X-ray fluxes. 

On the bases of these analyses, we assume that the 



WD is embedded in very hot, highly ionized gas that 
scatters supersoft X-ray photons emitted from the WD. 
In other words, the X-ray emitting region is not a point 
source but is substantially extended around the WD, as 
illustrated by the blue circular region in Figure [2] We 
also add circumbinary matter distributed as in Figure [2j 
in which the RG swelling over its Roche lobe has a long 
spiral tail of neutral hydrogen (see Section [4] for details). 

3. LIGHT CURVES OF V AND 7" VARIATIONS 

We assume a binary configuration as in Figure [T] to 
calculate the V and / light curves. As we already dis- 
cussed in Section |2~T| it is very likely that the RG almost 
fills its Roche lobe. We assume a circular orbit with an 
inclination angle of i = 90° and WD and RG masses of 
Mwd = 1.0-1.35 M Q and M RG = 0.8-2.0 M Q , respec- 
tively. We further assume the RG radius to be constant, 
although the RG is pulsating as a semiregular variable. 
This approximation is good if the amplitude of the ra- 
dius pulsation is small. In the symbiotic nova PU Vul, 
the companion RG is also a semiregular variable, and the 
amplitude of its rad ius pulsation is est imated to be 3-7% 
by eclipse analysis (|Kato et al. II2012T ). Because the am- 
plitude of V variation in SMC3 (~ 0.15 mag) is much 
smaller than that in PU Vul (~ 1 mag), a constant RG 
radius is a good approximation. 

The WD has hydrogen shell burning so its surface is 
hot enough to emit supersoft X-rays. The hemisphere 
of the RG toward the WD absorbs high-energy photons 
emitted from the WD, and this side is heated (see Fig- 
ure [TJ. The heated surface emits blackbody radiation at 
a local temperature calculated from the energy balance. 
In contrast, in the opposite hemisphere, the temperature 
remains low at the original surface temperature, Trg, 
which we take as a parameter. We neglected the tem- 
perature change due to RG pulsation. The RG surface is 
divided into small patches (N — N r x Ng = 64 x 32) and 
the energy balance is calculated at each patch to obtain 
its local temperature. Then, we sum them up to obtain 
the total V and I magnitudes of the RG at a specified 
orbital phase. A d etailed description of ou r calculation 
method is given in IHachisu fc Katol (l200llh 

The WD and its hot thin nebula (blue region in Figure 
[2) are extremely hot (T~5x 10 5 K: IQrio et alll2007h . 
so their emission contributes little to the V and I bands. 
We considered the hot nebula surrounding the WD to be 
optically thin (see the discussion in Section [5]) , so it is 
transparent in the / and V bands. We also assumed that 
the spiral tail is neutral and transparent in these I and 
V bands, as described in Section POl 

However, the cool tail may contribute to the I mag- 
nitude. As the standard model in this work, we deter- 
mined the original RG temperature (i.e., that of the non- 
irradiated surface) by fitting our V model light curves 
with the observed V magnitude. The resulting / light 
curve is darker than the observed one and does not fit 
well, as seen below. Therefore, we added a nebular con- 
tribution in the I band in order to fit our model light 
curve with the observed I band light curve. We suppose 
that this nebular contribution is due mainly to the cool 
RG tail. 

We assumed a steady hydrogen-shell-burning WD of 
bolometric luminosity Lwd = 3.77 x 10 4 Lq for Mwd = 
1.35 M Q , L wd = 3.59 x 10 4 L Q for M W d = 1.2 M Q , 
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HJD-2450000 

Fig. 3. — I and V light curves of SMC 3. Solid curves denote our 
model light curves for a 1.35 M WD and 0.8 M RG. The original 
temperature of the RG (i.e., the temperature on the non-irradiated 
surface), Trq, is the only parameter determined from light curve 
fitting in V band. We add a nebular emission of mj = 14.3 to 
fit the / band data. The / and V data were taken from OGLE 
II /III (red open circles) and MACHO (small green open circles). 
Our light curve model reproduces well the long-term variation in 
the V magnitude as well as the very small dynamic range (~ 0.1 
mag) of / magnitude variation. Irregular deviation from our model, 
especially in the I magnitude, may originate from the RG pulsa- 
tion, possibly in a higher overtone mode, which our model does not 
include. 

Lwd = 3.34 x 10 4 L Q for M WD = 1.1 M , and L WD = 
3.02 x 10 4 Lq for M WD = 1.0 M . We also assumed the 
RG mass to be Mr G = 0.8, 1.0, 1.5 and 2.0 M s . SMC3 
is located in th e direction of a dense H I cloud of iVn = 
4 x 10 21 cm- 2 (|Kahabka fc Pietsch Ifl996h , but SMC3 is 
in front of this cloud, because the column density derived 
from spectral analysis of supersoft X-rays is much smaller 
than this value {N-a ~ s everal xlO 20 cm" 2 : lOrio et al.~l 
[2007HSturm et al. H201lh . Therefore, we adopted E(B - 
V) = 0.099, a distance modulus of 18.91 (a distance of 
60.6 kpc), M = 0.179 and Ay = 3.1 x E(B - V) = 0.31 
after IHilditch et al~l (|2005T ) . 

Figurels hows the OGLE-II (JUdalski et al. lfl997h and 
OGLE-Ilfl (|Udalski et al. Il2008h I and V bands as well 
as the MACHC0 blue band light curve. The MACHO 
light curve was shifted in magnitude to match OGLE's 
V magnitude. This figure shows long-term variations in 
the / and V magnitudes, both of which are superposed 
with short pulsational variations of the RG. 

Our model light curves are shown in Figure [3] for a 
1.35 M Q WD and a 0.8 M RG. When the RG is in 
front of the WD (<fi ~ 0.0), we see only the original sur- 
face of the RG, where the magnitudes are faintest. On 
the other hand, when the RG is on the far side of the 
orbit {(f> ~ 0.3 - 0.7), we see the irradiated bright hemi- 
sphere of the RG. Thus, the light curve shows a wide dip. 
Secondary minima exist, as clearly shown in our model 
/ light curve. They are due to ellipsoidal variation (a 
double peak in one orbital period) of a lobe-filling com- 
panion. Here we assume that only the RG contributes to 
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Fig. 4. — Same as Figure [3] but for a 1.0 A/q RG. We added 
nebular emission of mi = 14.4 to fit the / band data. 
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Fig. 5. — Same as Figure [3] but for a 1.5 Mq RG. We added 
nebular emission of mi = 14.7 to fit the / band data. 



the V band emission, not the nebulae, the contribution 
of which is discussed later. We found that the RG sur- 
face temperature of the non-irradiated side Tr,g = 3460 
K yields a good result in fitting with the V light curve. 
To fit the I band light curve, we had to assume a nebular 
contribution of m\ = 14.3. 

If we assume no nebular contribution to /, we need to 
assume a higher temperature, Trg = 3700 K, to obtain 
a good result. This may be explained if we are observing 
different "photospheres" of the RG envelope at different 
wavelengths. When we observe the I magnitude, we may 
see a deeper region where the temperature is higher. An- 
other possible explanation is contamination by emission 
from the cool tail. The tail may contribute to the I band, 
but not the V band, reducing the temperature difference 
between different bands. Although the / band is contam- 
inated by semiregular pulsations, our / light curve ade- 
quately reproduces the very small dynamic range (~ 0.1 
mag) of the observed light curve. 
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Fig. 6.— Same as Figure |3] but for a 2.0 M RG. We added 
nebular emission of mj = 15.0 to fit the / band data. 

Figures H]-|6] show similar model light curves for a 1.35 
M Q WD with different RG masses, i.e., M RG = 1.0, 1.5, 
and 2.0 M Q . For a higher RG mass, we need to assume a 
lower RG surface temperature in order to obtain a sim- 
ilar magnitude. For a more massive RG, we obtained 
a larger separation a, and thus a larger effective Roche 
lobe radius (see Table [lj. Therefore, we had to assume 
a lower temperature in order to obtain a similar magni- 
tude. For each model, we assume a nebular contribution 
of tot=14.3-15 mag. 

As in the previous section, the RG pulsation is not a 
fundamental mode, but possibly exists in higher over- 
tones, which is consistent with the irregular pulsation 
behavior at small amplitude (~ 0.05 mag). These pul- 
sations are not included in our theory, so we cannot re- 
produce the complicated I band variations. Instead, it 
is important that our light curve model reproduces the 
upper and lower limits of the / magnitude around the 
mean value by the same model that reproduces the long- 
term variation in the V band. In this sense, we may 
conclude that our model is consistent with both the V 
and / observations. 

In Figures [3] [6l the V magnitude is reproduced well by 
all of our theoretical models, whereas the / magnitude 
may be better reproduced by a low-mass RG companion. 
As shown in the previous subsection, Trg = 3700 ± 100 
K is quite consistent with the pulsation theory, so we 
prefer a 0.8 M RG model to the others. 

We also calculated V and / light curves for different 
masses of the WD, i.e., Mwd = 1-2, 1.1, and 1.0 M Q , 
which are more or less similar in both the V and / mag- 
nitudes (Figures [3 and[9j respectively). Because both 
the I and V light curves are contaminated by irregular 
RG pulsations, we cannot determine which is the best 
fit to the observation. Thus, we could not determine the 
WD mass (or constrain the range of the WD mass) using 
only these light curve fittings. 

The outer part of the hot thin nebula around the WD 
could contribute to the V band because it has a large vol- 
ume even though it has very low emissivity. Part of the 
cool RG tail could also contribute to the V band if it is 
ionized by irradiation. Thus, we calculated an additional 
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Fig. 7.— Same as Figure [3] but for a 1.2 A/ WD. 
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nebula = 1 4.3 mag 
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Fig. 8.— Same as Figure [3] but for a 1.1 A/ WD. 
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nebula = 1 4.3 mag 
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Fig. 9.— Same as Figure [3] but for a 1.0 M Q WD. 
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M V1D =1.35M S M RG =0.BM s T RG =3300 K P„ rb =1634 days (m-M) v =19.2 
nebula=16.5 mag (30%) 
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Fig. 10. — Same as Figure [3] but with nebular emission of 
my = 16.5 and mj = 14.0 (thick solid blue lines). Thin solid ma- 
genta line denotes the original V magnitude (without nebulosity) 
shifted upward by —0.4 mag. The nebular emission of my = 16.5 
corresponds to about 30% of the original luminosity. The calcu- 
lated eclipses become too shallow to reproduce the observation. 

light curve model in which the nebular contribution to 
the V band is as much as 30% of that of the RG. Figure 
[10] shows such a case for a pair of 1.35 Mq WD and 0.8 
Mq RG. For example, if the nebular emission amounts 
to 30% of the RG luminosity in the V band, we obtain 
a lower temperature of Trg = 3300 K for the RG com- 
panion. The total V magnitude is shown by a blue line. 
To fit our model with the / magnitude at the same RG 
temperature, we must increase the nebular contribution 
from mi = 14.3 to mi = 14.0. If we compare this model 
with Figure [31 the amplitude between minima and max- 
ima becomes smaller in the V band, which is inconsistent 
with the observations. Therefore, we may say that the 
nebular contribution is probably less than my > 16.5. 

In many symbiotic stars, the / and V band magni- 
tudes are strongly contaminated by nebular emission 
(iNussbaumer fc Vogel Hl987tlSkopal et al. lll998tlSkopal I 
I2005t iSkopal et al. Il2009i rin SMC3, however, the com- 
panion RG is a very bright AGB star. This explains why 
the nebular contributions in / and V are rather small. 

In our V light curve fitting, we use the ephemeris 

MJD min ,v = 49280 + N x (1634) days. (2) 

This is consistent wi th those of JD min = 2449360 + N x 
(1600 ± 140) days bv lKahabkal (I20041). and MJD n ^ B = 
(49242±9)+iVx (1647±24) davs bv lSturm et al. I (|201lh . 

Periodic variations similar to those of SMC3 were ob- 
served in YY Her, a symbiotic binary of P rb = 590 days. 
The RG fills its Roche lobe, and it s hemisphere toward 
the W D is heated by irradiation. iMikolaicw ska et al. I 
(2002) showed that the periodic change in the UBVRI 
magnitudes can be described by a combination of ellip- 
soidal and sinusoidal changes. The RG radiation domi- 
nates in longer wavelength bands among U BVRI, so the 
variation can be represented well by ellipsoidal changes. 
On the other hand, the nebular continuum dominates 
in shorter wavelength bands, so the light curve shows 
a single minimum in one orbital period. These charac- 
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Fig. 11. — (a) X- ray light curve of SMC3. Data are taken from 
ISturm et af~l feoTlT ). Filled triangles: ROSAT PSPC. Open trian- 
gles: ROSAT HRI. Open circles: Chandra. Open squares: XMM- 
Newton. Stars: Swift. The supersoft X-ray flux (0.2 - 1.0 keV) is 
in units of 10 — 12 erg cm -2 s — 1 (left ordinate). Solid lines indicate 
model light curves of X-ray eclipses; Model 1 for -Rrg = 245 Rq 
(thick red line), Model 2 for R RG = 199 Rq (thin black line). 
Right ordinate indicates the ratio of the calculated X-ray flux to 
the original value due to partial/total eclipse by the RG and its 
tail. Vertical dashed line indicates the central time of X-ray total 
eclipses for both Model 1 and Model 2 (0 = 0.075). (b) Height 
of the RG and its tail from the orbital plane. Model 1 (thick red 
line). Model 2 (thin black line). Red and black dotted lines denote 
the RG surface in Models 1 and 2, respectively. Horizontal dashed 
line indicates the radius of the X-ray-emitting sphere around the 
WD. 

teristic properties are already shown in our model light 
curves for SMC3. Figures [SHU show that the I magnitude 
variations have double peaks whereas the V light curves 
display a single eclipse because the secondary minimum 
is compensated by extra radiation from the irradiated 
side of the RG. 

4. X-RAY LIGHT CURVE 

4.1. A Wide X-ray Eclipse 

Figure [TTk shows the X-ray fluxes of SMC3 against the 
orbital phase superposed for five orbital period s, the data 
for which are taken from ISturm et aT~l (|201lD . assuming 
the ephemeris and period determined by Equation ([2]) 
in Section [3] The same data are repeated for two or- 
bital cycles. This X-ray light curve shows the following 
characteristic properties. 

(1) : The X-ray eclipse is so wide that it spans A</> ~ 

0. 15-0.3 of the orbital phase represented by Equa- 
tion ©. 

(2) : The central time of X-ray eclipses is at (f> ~ 0.1, 

1. e., ~ 150 days after <f) = 0.0 as determined by 
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Equation @, i.e., as determined from the long- 
term V variation. 

(3): The eclipse has an asymmetric shape; the ingress is 
much steeper than the egress. 

Regarding property (1), the eclipse duration is not 
unusually long among symbiotic stars. The dura- 
tion depends on the wavelength because of absorp- 
tion/scattering by surrounding material (RG winds) out- 
side the RG photosphere. For example, in the symbi- 
otic nova PU Vul (P rb = 13.5 years), optical eclipse 
lasts A(f) — 0.07, which is longer than the RG ge- 
ometrical eclipse A. 6 = 0.044 because of TiO atmo- 
spheric absorption (|Kato et al. 1 12012( 1. At UV wave- 
lengths region, its eclipse duration is even longer be- 
cause of Rayleigh scattering as Acj> > 0.14 (up to 0.27) 
(|Tatarnikova fc Tatarnikov H2009I; iKato et al. 1 1201 111 . 

Regarding property (2). iSturm et al. I (|2011 1 already 
reported that the X-ray minimum comes 140 days af- 
ter the V magnitude minimum. They noted this as B 
instead of V, but the MACHO b lue is closer to yellow 
than B (jBessel fc Germany 111999ft . They attributed this 
difference between the minima to a possible error in the 
V magnitude minimum due to the short observing period 
of MACHO (about 1.5 orbital periods). We see, however, 
a clear difference between the X-ray mid-eclipse and V 
magnitude minima because we determine the ephemeris 
from 3.5 orbital periods as shown in Figures [3HSl 

All of these properties are consistent with an eclipse by 
an M giant companion swelling over the Roche lobe on 
the trailing side of the binary orbit, as wa s already been 
suggested in the symbioti c star SY Mus (jDumm et al. I 
[199a iPereira et al. Ifl995h . The RG cool winds form a 
thick envelope, which tends to flow out mainly backward 
(on the trailing side of the RG) owing to orbital mo- 
tion. As a result, SY Mus shows a wide eclipse in the 
UV continuum (</> ~ 0.9 to ~ 1.18, i.e., A<j) = 0.28), 
which is much wider than the geometric eclipse by the 
RG companion itself (Ac/> = 0.074). Moreover, the UV 
1325 A continuum shows an asymmetric eclipse in which 
the ingress is steeper and closer to the mid-eclipse than 
the egress. This property is th e same as property (3) 
of SMC3. iDumm et al~l fjl999T ) analyzed the UV spec- 
trum and concluded that these eclipse properties are nat- 
urally explained by obscuration of the UV flux by a thick 
neutral hydrogen cloud (Nn from 10 23 cm~ 2 to > 10 25 
cm~ 2 ) surrounding the RG. 

Because the X-ray eclipse of SMC3 shows similar prop- 
erties, we suppose that it is also caused not only by the 
RG companion itself, but also by a thick atmosphere 
surrounding the RG. In other words, properties (l)-(3) 
above are evidence of eclipses caused by a companion 
star swelling over the Roche lobe on the trailing side of 
the binary orbit (see Figure [5]) . 

4.2. Circumbinary Configuration 

For the X-ray emitting region, we simply assume a 
spherical hot nebula of radius i?NE centered on the WD. 
The WD itself emits stronger X-rays than the surround- 
ing nebula, but the WD is believed to be always obscured 
by a neutral hydrogen cloud surrounding the binary in 
the orbital plane. Therefore, we observe only the ex- 
tended X-ray source (blue region in Figure [2]) which is 



always partially occulted by the trailing tail from the 
RG. We will discuss the X-ray emitting region in Section 

El 

We assume that the surface brightness of the X-ray 
emitting region is constant. This may not be a good 
approximation, especially if the hot nebular gas origi- 
nates in the wind from the WD, the density of which 
decreases as r~ 2 , where r is the distance from the WD. 
lOrio et al~l (|2007l) , however, suggested that their spec- 
tral analysis with CLOUDY favors a constant density 
distribution over an r~ 2 law. In that case, the surface 
emissivity (emissivity per unit area of the hot nebula in 
the line of sight) is close to constant rather than decreas- 
ing with increasing power of the distance from the WD. 
We suppose that if the X-ray flux is due to Thomson 
scattering with an optical depth of about 1 rather than 
thermal emission, the surface emissivity may not depend 
strongly on the distance from the WD. In this paper we 
do not discuss whether the X-ray flux is due to Thom- 
son scattering, thermal emission of the hot plasma itself, 
or another cause, but, simply assume a uniform surface 
emissivity (a circular disk with a constant emissivity) as 
a first approximation. 

Many numerical calculations show that cool RG 
winds (and also hot WD winds if present) form 
spira l structures around t he b in ary dSawada et al~1 
1981 IGawryszczak et al~l 12001 Kim fc Kim I 1 20071: 
Edgar et al. I 120081: IMohamed & Podsiadlowski I bold 
and references therein). Ejected matter is highly aspher- 
ical, ha ying a tendency to concentrate toward the orbital 
plane (iBisikalo et al 1 120001: IGawryszczak et al~l l2002t 
lEdgar et al.T l2008: M ohamed fc Pod siadlow ski 1 12012ft . 
In optical spectra of SM C3 ((Morgan H992tlMiirset et all 
119961: lOrioet al. 112007ft . Raman-scattered O VI lines in- 
dicat e a significant amount of neutral material ((Schmid I 
11996( 1. Therefore, we suppose that SMC3 has a tail of 
neutral material concentrated in the orbital plane. 

Figure [5] shows a schematic picture of our circumbinary 
model of X-ray eclipses. The RG is swollen to almost 
twice in its size in the direction opposite to the orbital 
motion and is followed by a long, geometrically thick tail 
surrounding the binary in the orbital plane. The shape 
of the RG and its tail is drawn schemat ically using in- 
forma tion fr om numerical calculations by iSawada et al. I 
(|1986ft and IMohamed fc Podsiadlowski I (|201j), which 
show narrow spiral tails flowing from both stars. Mo- 
hamed & Podsiadlowski's 3-D calculation showed that 
the outflowing matter is highly aspherical and concen- 
trated in the orbital plane. The details of the mass dis- 
tribution in the orbital plane do not affect our X-ray 
light curve model because we assume a total eclipse with 
an inclination of i = 90°. Therefore, only the heights 
(perpendicular to the orbital plane) of the RG and tail 
are important in our X-ray light curve modeling, which 
we use to reproduce the X-ray light curve. We suppose 
that this tail is composed of neutral hydrogen, which is 
opaque only to supersoft X-ray photons but transparent 
to optical photons. As shown later, the radius of the X- 
ray emitting region i?NE = 100 Rq, which is about half 
of the RG radius. This ratio is roughly consistent with 
the mass distribution aroun d the WD in the 3-D calcula- 
tion of binary mass outflow (Mohamed & Podsiadlowski 
12012ft . although their binary parameter differs from ours. 

Such neutral hydrogen material, opaque to X-rays but 
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transparent to optic al photons, was r eported in the re- 
current nova U Sco (|Ness et al. II2012D . In the supersoft 
X-ray phase followed by the 2010 outburst, U Sco tem- 
porarily showed irregular rectangular-shaped variations 
in the X-ray flux, whereas the UV and blue l ight curves 
showed a smooth eclipse. iNess et aT~l (j2012T ) explained 
this X-ray variation as occultations of an extended X- 
ray region due to Thomson scattering by several blobs 
of neutral hydrogen moving in front of the X-ray emit- 
ting region. These blobs obscured supersoft X-rays but 
were transparent to UV/blue optical photons. In their 
analysis, the X-ray emitting region extends beyond the 
binary orbit (the orbital period is P or b = 1-23 days, and 
the separation is a = 6.5i? Q ), and the blob sizes (0.73- 
1.1 Rq) are comparable to that of the companion MS 
star (i?MS = 2.1 Rq). Because SMC3 is a wide binary 
(P orb = 1634 days, a = 700-1000 R & as tabulated in Ta- 
ble [IJ and its SSS phase lasts much longer than that of U 
Sco (which lasts only 15 days), the circumstellar matter 
could be much more massive and extend to a much larger 
size than in U Sco. 

In this paper, we assumed that the eclipse is essentially 
total; i.e., the X-ray emitting region (radius -Rne) is com- 
pletely covered by the RG, although weak X-ray fluxes 
< 1/20 of the X-r ay high state) have been observed 
Sturm et al. 1120111) . We regard these weak X-rays as 
light scattered by tenuous plasma widely distributed be- 
yond the radius -Rne (= 100 Rq). This assumption does 
not affect our light curve modeling because the scattered 
light is very faint. 

4.3. Model X-Ray Light Curve 

Our model, which reproduces the observed X-ray vari- 
ation, is shown in Figure [TT] The orbital phase (f> = 0.0 
is defined as when the RG is in front of the WD. The 
theoretical X-ray light curve is calculated from the area 
of the X-ray emitting region occulted by both the RG 
and its tail at a given orbital phase. Model 1 (denoted 
by a thick red line) assumes a greatly swollen RG with 
a thick (i.e., taller) tail and reproduces the lowest X-ray 
flux ever observed. On the other hand, Model 2 (denoted 
by a thin black line) assumes a smaller RG with a thin 
tail and reproduces the highest ever observed X-ray flux. 
The geometrical thickness of the tail is shown in Figure 
II lb . and its minimum height is listed in Table [2] 

We found that an RG radius of 200-280 Rq is consis- 
tent with the ingress steepness, the value of which de- 
pends weakly on -Rne and the height of the tail. The 
nebula size of -Rne = 100 Rq may be a bit large, but it 
cannot be a point source compared with the RG because 
the ingress shows a finite slope that is not perpendicular. 

Outside of total eclipses (i.e., at the orbital phase 
4> ~ 0.3-0.8), the X-ray data are scattered, as shown 
in Figure II lb : this represents real variations between 
different orbital cycles or is attributed to possible er- 
rors in the conversi on factor between different satellites 
(ISturm et al. ll20Ol . If this is a real variation, different 
X-ray fluxes are attributed to different shapes (heights) 
of the RG tail. Because the RG is pulsating with a period 
of 110 days, the mass- loss rate of the cool wind could vary 
with time, in addition to exhibiting possible long-term 
irregular variations. Thus, the RG tail naturally has dif- 
ferent thicknesses (heights) at different times. (Note that 
periodic variations in the tail shape are observed in a 3-D 
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Fig. 12. — Photospheric radius versus photospheric temperature 
of our envelope models for WDs with masses of 0.6, 0.8, 1.0, 1.25, 
and 1.35 Mq. The chemical composition of the envelopes is as- 
sumed to be X = 0.7, Y = 0.298, and Z = 0.002. Red filled circles 
indicate places at logg = 7.0, 7.5, 8.0, 8.5, and 9.0. Optically 
thick winds occur in the region left of open square. No optically 
thick winds are accelerated by the 0.6 and 0.8 Mq WDs because 
of low metallicity in the SMC. Short horizontal bars indicate the 
Chandrasckhar radii corresponding to each WD mass. 

calculation of the mass outflow from Mira with a peri- 
odica lly changing mass-loss rate having a period of 332 
days: iMohamed fcP odsiadlowskil 120121) . This explains 
the short-term variation in the X-ray flux as well as the 
differences among orbital periods. Note that the combi- 
nation of the RG size and geometrical thickness (height) 
of the tail is not unique. We determined the tail height 
relative to -Rne = 100 Rq. If we adopt a smaller Rne, for 
example, a similar X-ray light curve is obtained with the 
tail height multiplied by -Rne/(100 Rq) to the present 
models. The detailed distribution of circumbinary mat- 
ter in the azimuthal and radial directions on the orbital 
plane does not affect our X-ray eclipse light curve model 
because only the height of the tail at a given orbital phase 
is important for determining the X-ray flux. 

Figure [TTk also indicates the central time of the X-ray 
total eclipse by a dashed line, according to the ephemeris 
determined from the V magnitude minimum. The X-ray 
mid-eclipse comes 123 days after (f> — 0.0 determined 
from the minima in the V light curve. 

5. DISCUSSION 

Figure Q2] gives the photospheric radius versus the pho- 
tospheric temperature for WDs of different masses with 
steady hydrogen-shell-burning. The chemical composi- 
tion of the envelope is assumed to be X — 0.7, Y = 
0.298, and Z = 0.002 Th e nu merical method is th e 
same as that in lKato I ([1999) and lNomoto et al. I (|2007f) . 
In the region to the left of the open squares, optically 
thick winds occur and absorb supersoft X-rays. In the 
region to the right of the open squares, the WD envelope 
is in hydrostatic equilibrium, which corresponds to a SSS 
phase. The lowest end of each line (i.e., the smallest WD 
radius) corresponds to the point at which hydrogen burn- 
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ingjs extinguished. lOrio et al~l (|2007l) and lSturm et al. I 
(|201lD obtained an X-ray temperature for SMC3 as soft 
as 30-50 eV. If this temperature corresponds directly to 
the photospheric temperature, the radius of the photo- 
sphere is less than 0.1 Rq, much smaller than the radius 
of the hot nebula, i?NE = 100 Rq. This is consistent 
with our assumption that the WD is always totally ob- 
scured by the RG and its long tail, and we see only X-rays 
scattered by the hot nebula. 

This figure also shows the positions of the WD en- 
velopes corresponding to surface gravities of logg=7.0, 
7.5, 8.0, 8.5, and 9.0, where g = GM WD /Rl h is the 
photospheric gravity. There are no envelope solutions 
corresponding to log t/=8.5 and 9.0 for Mwd < 1-0 Mq. 
IQrio et al~l (|2007[ ) suggested a large gravity of log g ^ 9.0 
and a law temperature of 36-47 eV on the basis of X-ray 
spectral fitting. Such a large gravity corresponds, in our 
WD envelope models, to a very massive WD (> 1.3 M©) 
with a small radius, and also a much higher temperature 
(> 100 eV) as shown in Figure [T^l The origin of this 
discrepancy is unc lear. 

On the basis of IQrio et al"~l (|2007D . we assumed that 
the WD is surrounded by a very hot nebula, which scat- 
ters X-rays from the WD. We estimate the density of the 
hot plasma from the coronal lines. In principle, the elec- 
tron density in the coronal region should be close to the 
critical density; e.g., the critical den sity for [Fe X] 6374 is 
4.8 x 10 9 cm" 3 (jNaeao et al. |[200l . In SMC3, both the 
[Fe X] 6374 and [Fe IX] 7892 lines are relatively strong 
and thus trace regions of high electron density, N e ~ 10 9 - 
10 10 cm -3 . For a spherical emitting region of radius 
Rne = 100 Rq (fully ionized constant-density gas) the 
total mass of the nebula is estimated to be roughly 10 -9 - 
10 -8 Mq. The optical depth of the region in the line of 
sight is r ~ npr — 0.004-0.04 for the electron scattering 
opacity. Therefore, this region is optically thin, which is 
consistent with our assumption. However, these values 
are highly uncertain because of the simplified assump- 
tion. 

We also estimate the temperature of this region. These 
coronal lines, [Fe X] and [Fe IX], have high ionization 
potentials (IPs), of 234 eV and 262 eV, respectively, 
and their presence indicates that the ionized matter is 
at a high temperature. Using the formula T/(10 3 K) 
~ ma ximum observed IP (eV) (|Miirset fe Nussbaumer I 
[1991 . we obtain a temperature for the ionizing source of 
T =(2-3) xlO 5 K. This plasma temperature is consistent 
with the WD temperature estimated from X-ray spec- 
trum analysis, 36-47 eV (4.2-5.5 xlO 5 K) (jOrio et al. I 
[20071) . 

6. CONCLUSIONS 

We estimated the radius of the RG and confirmed that 
it almost fills its Roche lobe. On the basis of this pic- 
ture, we calculated the variations in orbital brightness 
including an irradiation effect and the geometrical ef- 
fect of a tidally distorted companion. Our model reason- 
ably reproduces the observed long-term / and V light 



curves of SMC3. This again confirms that the RG al- 
most fills its Roche lobe. Thus, the mass transfer rate 
from the RG to the WD could be high enough (> a few 
times 10 -7 Mq yr _1 ) to ma intain steady hydrogen -shell- 
burning on the WD (e.g., iNomoto et al. I [20071 ). even 
though SMC3 is a long orbital period binary (P or b = 4.5 
years). This is quite consistent with the fact that SMC3 
has been a steady SSS for at least 20 years. 
Our main results are summarized as follows. 

1. We presented I and V light curve models of SMC3. 
The observed magnitude variations can be repro- 
duced by an irradiation of the RG by the hot WD 
as well as geometrical effects of the tidally dis- 
torted RG. This confirms that the RG component 
almost fills its Roche lobe, suggesting high mass 
transfer rates onto the WD, which maintain steady 
hydrogen-shell-burning on the WD. 

2. We also presented an X-ray eclipse model in which 
an X-ray-emitting region is eclipsed by an RG com- 
panion swelling over the Roche lobe on the trailing 
side of the binary orbit. The X-ray-emitting re- 
gion is highly ionized and extends to as large as 
~ 100 Rq, scattering supersoft X-rays from the hot 
WD. This configuration reasonably reproduces the 
asymmetry of the X-ray eclipses, that is, steeper 
ingress than egress and the occurrence of the cen- 
ter of the X-ray eclipse about 120 days after the V 
minimum, as well as a very wide X-ray eclipse. 

3. We also showed that our circumbinary matter con- 
figuration is consistent with various observed char- 
acteristics; our configuration is composed of a very 
hot nebula surrounding the WD and circumbinary 
matter (a spiral tail) of neutral hydrogen that al- 
ways totally obscures the WD but partially occults 
the extended X-ray emission region. These con- 
figurations may be the key to the weakness of the 
supersoft X-ray flux in symbiotic systems. 
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